Skip to content

Conversation

@njeffery
Copy link
Contributor

@njeffery njeffery commented Mar 28, 2025

PR checklist

  • Short (1 sentence) summary of your PR:
    Corrects bug introduced in Icepack v1.5.0 that causes negative silicate values.
    See: BGC bug in Icepack v1.5.0 E3SM-Project/E3SM#7127

  • Developer(s):
    @njeffery

  • Suggest PR reviewers from list in the column to the right.

  • Please copy the PR test results link or provide a summary of testing completed below.
    Tested in a fully-coupled biogeochemistry simulation. Simulation page and comparisons with a control run are here:
    https://acme-climate.atlassian.net/wiki/spaces/HESF/pages/5121212417/20250314.v3.LR.CBGC.FrazilEcosys.noSpinup.pm-cpu

  • How much do the PR code changes differ from the unmodified code?

    • bit for bit (for physics only runs)
    • different at roundoff level
    • more substantial for sea ice bgc
  • Does this PR create or have dependencies on CICE or any other models?

    • Yes
    • No
  • Does this PR add any new test cases?

    • Yes
    • No
  • Is the documentation being updated? ("Documentation" includes information on the wiki or in the .rst files from doc/source/, which are used to create the online technical docs at https://readthedocs.org/projects/cice-consortium-cice/.)

    • Yes
    • No, does the documentation need to be updated at a later time?
      • Yes
      • No
  • Please document the changes in detail, including why the changes are made. This will become part of the PR commit log.
    The ice.log files with Icepack v1.5.0 with active BGC has warning messages after about 4 years indicating negative silicate. This is because silicate limitation is not implemented correctly. The fix removes an incorrect if block around the silicate limitation functions, which was mistakenly added in Icepack v1.5.0. With this fix, no warning messages were observed over the simulation (~36 years).

@eclare108213 edited to move template from comment below to top of PR.

eclare108213 and others added 2 commits November 21, 2024 10:42
Merge CICE-Consortium/Icepack main (v1.5.0) to E3SM-Project/Icepack main
Fixes a bug introduced in V1.5.0.
Silicate limitation should not be wrapped in an if statement for iron tracers.

nonBFB with active BGC.  BFB for all else.
@njeffery njeffery marked this pull request as draft March 28, 2025 19:45
@njeffery
Copy link
Contributor Author

@apcraig: I made a PR to E3SM-Project/Icepack identical to this one -- E3SM-Project#39. Let me knowif you'd like me to submit the PR here, keep as draft, or close.

@eclare108213
Copy link
Contributor

@apcraig: I made a PR to E3SM-Project/Icepack identical to this one -- E3SM-Project#39. Let me knowif you'd like me to submit the PR here, keep as draft, or close.

Since v1.5 there are only a few updates to Icepack. This is the only one that would affect E3SM -- the others are machine ports and similar things not related to the column physics, plus one other bug fix for the case calc_Tsfc = .false.. My preference would be to merge this BGC bug fix here then pull the whole batch since v1.5 over to E3SM, to keep them in sync. But we can do them in parallel if there are concerns about doing it this way. @proteanplanet

@apcraig
Copy link
Contributor

apcraig commented Apr 1, 2025

Happy to merge this. I like the idea of pull the latest Icepack into E3SM, but that's an E3SM process issue. Looks like the PR template hasn't been filled out for this PR, that's something that should probably happen.

@njeffery
Copy link
Contributor Author

njeffery commented Apr 1, 2025

PR checklist

  • Short (1 sentence) summary of your PR:
    Corrects bug introduced in Icepack v1.5.0 that causes negative silicate values.

  • Developer(s):
    @njeffery

  • Suggest PR reviewers from list in the column to the right.
    @eclare108213 and @apcraig

  • Please copy the PR test results link or provide a summary of testing completed below.
    Tested in a fully-coupled biogeochemistry simulation. Simulation page and comparisons with a control run are here:
    https://acme-climate.atlassian.net/wiki/spaces/HESF/pages/5121212417/20250314.v3.LR.CBGC.FrazilEcosys.noSpinup.pm-cpu

  • How much do the PR code changes differ from the unmodified code?

    • [x ] bit for bit (for physics only runs)
    • different at roundoff level
    • [x ] more substantial for sea ice bgc
  • Does this PR create or have dependencies on CICE or any other models?

    • Yes
    • [ x] No
  • Does this PR add any new test cases?

    • Yes
    • [x ] No
  • Is the documentation being updated? ("Documentation" includes information on the wiki or in the .rst files from doc/source/, which are used to create the online technical docs at https://readthedocs.org/projects/cice-consortium-cice/.)

    • Yes
    • [x ] No, does the documentation need to be updated at a later time?
      • Yes
      • [x ] No
  • Please document the changes in detail, including why the changes are made. This will become part of the PR commit log.
    The fix removes an incorrect if block around the silicate limitation functions which was mistakenly added in Icepack v1.5.0.

@njeffery njeffery marked this pull request as ready for review April 1, 2025 22:25
@eclare108213 eclare108213 self-requested a review April 2, 2025 15:50
@apcraig
Copy link
Contributor

apcraig commented Apr 2, 2025

@apcraig apcraig merged commit b3b5127 into CICE-Consortium:main Apr 2, 2025
2 checks passed
Sign up for free to join this conversation on GitHub. Already have an account? Sign in to comment

Labels

None yet

Projects

None yet

Development

Successfully merging this pull request may close these issues.

3 participants